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Abstract 

The trajectories of differently shaped nanoparticles manipulated by atomic force microscopy are related to the scan path of the 
probing tip. The direction of motion of the nanoparticles is essentially fixed by the distance b between consecutive scan lines. Well- 
defined formulas are obtained in the case of rigid nanospheres and nanowires. Numeric results are provided for symmetric nano- 
stars. As a result, orienting the fast scan direction perpendicular to the desired direction of motion and reducing b well below the 
linear size of the particles turns out to be an efficient way to control the nanomanipulation process. 



Introduction 

Quite soon after its invention, it became clear that atomic force 
microscopy (AFM) could be used not only for maging but also 
for manipulating nano-objects [1,2]. This possibility has 
produced spectacular results and last, but not least, it has 
allowed the controlled manipulation of metal clusters on insu- 
lating surfaces [3] and even single atoms on semiconductors [4]. 
However, AFM manipulation tends to be time-consuming. A 
major issue is that nanoparticles are usually moved individually 
so that the AFM tip has to be properly positioned with respect 
to the particle every time. The tip is either placed on the side or 
on the top of the particle. Then the tip-particle interaction is 



increased (by varying the tip-particle distance or the amplitude 
of the tip oscillations) until the particle is detached from the 
substrate and moved in a direction which is determined by 
several factors such as the scan pattern, the surface structure 
and the geometry of both tip and particle. Predicting the direc- 
tion of motion of nanoparticles is very important, especially if it 
is desired to manipulate several particles at the same time. Here, 
we show that this is possible in simple cases of practical 
interest. Specifically, we assume that the AFM is operated in 
tapping mode (although some conclusions may be extended to 
contact mode), the particles are rigid and the frictional forces 
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between particles and substrate can be neglected when the parti- 
cles collide with the tip, but they are high enough to stop the 
particles immediately once contact with the tip is lost. The 
concentration of nanoparticles on the substrate is also supposed 
to be low enough to prevent multiple collisions in the manipula- 
tion. 

After a brief review of previous results on the manipulation of 
rigid nanorods, including nanospheres and thin nanowires as 
limit cases, we discuss symmetric nanostars as a prototype of 
more complex shaped particles. We show that in any case the 
angle of motion of the nanoparticles is precisely related to the 
distance b between consecutive scan lines. When the parameter 
b is sufficiently small, the particle tends to move perpendicu- 
larly to the scan direction. The exact relation between the angle 
of motion 6 and the parameter b depends on the particle shape 
and can be, in principle, determined analytically. Curiously, this 
has a certain analogy to the scattering of sub-atomic particles, 
whose angle of deflection 9 depends on the form of the scat- 
tering field and on the impact parameter b (i.e., the distance at 
which the particle would pass the center of the field in the 
absence of any interaction) [5]. 

Results 
The model 

We first consider a planar island whose profile is described by 
the function r = r(cp) in polar coor-dinates or, equivalently, by a 
multi-value function y = y(x) in cartesian coordinates. Assuming 
that the tip follows a raster scan pattern, the y coordinate of the 
tip varies as Yq = Nb, where TV is the number of the scan line 
and b is the distance between consecutive scan lines 
(Figure la). 




Figure 1 : (a) A sharp nanotip follows a raster scan pattern with 
consecutive scan lines separated by a distance b. The tip collides with 
a nanoparticle (here represented by a star-shaped island) at the loca- 
tion P. (b) In tapping mode the tip oscillates in the direction z perpen- 
dicular to the plane of the figure and applies an impulsive force F 
perpendicular to the island profile, (c) In contact mode the force F is 
directed along the x axis and the total force acting on the particle will 
be oriented as in tapping mode only if the static friction f can balance 
the component of F along the island profile. 



The island has a mass M and a moment of inertia / with respect 
to the normal axis z through its center of mass (COM). Here, we 
assume that the linear size of the island is much larger than the 
tip radius, so that the force F applied by the tip is concentrated 
at the point of contact P. We also assume that the island cannot 
be deformed or broken during the manipulation. In such a case, 
the position R = (X,Y) of the COM and the angle of rotation O 
of the island about the normal axis z evolve according to the 
equations of motion of a rigid body: 



F = M 



d_R 

dt 2 



(1) 



and 



d 2 <t> 

P dt 2 



(2) 



where rp defines the position of the point of contact P with 
respect to the COM. 

The direction of the force F depends on the operating mode of 
the AFM. In tapping mode the tip oscillates in the z direction 
with a frequency in the order of 100 kHz with an amplitude of 
some tens of nm. This corresponds to an average speed of some 
mm/s, which is well above typical scan velocities in AFM 
(normally in the order of 1 jim/s). Thus, the tip hits the particle 
almost vertically and the vector F is oriented perpendicularly to 
the island profile, i.e., at an angle a = p + 90° with respect to the 
x axis, where 



tanp: 



dy 
dx 



r' p sin (p p + r p cos (p p 
r' p cos (p p + r p sin cp p 



and r' is the first derivative of r(cp) with respect to cp (Figure lb). 
In contact mode the tip hits the particle along the x direction and 
the force F can be oriented as in tapping mode only if the static 
friction force f between tip and particle is high enough to 
prevent sliding along the island profile (Figure lc). 

Assuming that friction between island and substrate is also high 
enough to prevent any slippage of the island after a collision 
with the tip, Equation 1 and Equation 2 can be averaged over 
the short collision time At (in the order of 1/ /, with /~ 10 5 Hz 
being the oscillation frequency of the tip). This leads to the 
equations 



dY 
dX 



tan a , 



(3) 
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for the translation of the island and 

dO _Mr v sin(a-(p p ) 
dR I 



(4) 



for the rotation, which can be finally integrated over the total 
time of interaction between tip and particle (along the given 
scan line). If the nanoparticle is not flat, it is easy to see that the 
previous analysis is still applicable provided that the particle 
does not roll and that its shape is not cylindrical. 

Translation and wobbling of nanorods 

The manipulation of a rigid nanorod formed by a cylinder (with 
length L) and two hemispherical caps (with radius a) is particu- 
larly instructive. Here, any possible rolling can be ignored and 
we can distinguish between two types of collision: (a) The tip 
touches the cylindrical core of the nanorod ("core" collision), 
(b) The tip touches one of the two hemispherical ends of the rod 
("cap" collision). In case (a) the equations of motion of the 
nanorod can be written in the form [6] 

dO 2M(Y 0 -Y-acosO) 
~dY~ /sin 20 ' 



and 



dY 
dX 



= - cot O 



In the case (b): 



dO 
dY 



ML sin(0 - a) 
21 sina 



(5) 



and 



dY 
dX 



= tan a . 



(6) 



In general, both core and cap collisions occur along each scan 
line and only numerical solutions are possible. However, a 
complete solution can be found in two important cases: The 
manipulation of a nanosphere of radius a (L = 0) and that of a 
thin nanowire of length L (a = 0), where only cap collisions or 
core collisions, respectively, occur. In the case of a nanosphere, 
Equation 5 and Equation 6 can be integrated leading to the 
following result [7]. The direction of motion of the sphere 
forms an angle 9 with respect to the x axis (fast scan direction) 
given by 



tan9 : 



a [cos a 0 + log tan ( a 0 /2 )] 



(7) 



The quantity ao is the impact angle between tip and sphere 
(with the exception of the very first collision) and is given by 



a n =arcsin 1- 



In the case of a nanowire, the average direction of motion is 
well-defined and is given by the sim-ple formula [6] 



cos 9 = — . (8) 
L 



The wire oscillates perpendicularly to this direction: 
0 = 9+90°. 

Thus, Equation 7 and Equation 8 show that the directions of 
motion of nanospheres and nanowires manipulated by AFM in 
tapping mode are completely determined by the distance b 
between consecutive scan lines or, equivalently, by the density 
of scan lines lib. The functions of Equation 7 and Equation 8 
are plotted in Figure 2. In both cases 9(Z?) decreases with 
increasing b until the particle is lost when b > a or b > L. 
Furthermore, the angle 9 — > 99° when b — > 9. Numerical simula- 
tions show that similar conclusions are also valid for arbitrarily 
thick nanorods [6], although simple analytic expressions cannot, 
in general, be derived. 




Figure 2: Angle of motion 9 of a nanosphere (solid curve) and a 
nanowire (dashed curve) as a function of the distance b between 
consecutive scan lines. The parameter b is expressed in units of the 
sphere radius a and wire length L respectively. 
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Star shaped islands: Rotational effects 

As a next step we extend our analysis to more complex shapes. 
We consider star-shaped islands, whose profile is described by 
the function 

r(cp) = a(l + sin 2 kcp) . 

The number of branches in the island is denoted by 2k. For 
instance, k = 3 in Figure la. It is inter-esting to observe that 
both the moment of inertia / and the area A of the island are 
independent of k: 

2n r(q>) 

I = ojd(p | r 3 dr = naa 4 , 

o o ^56 

A = \ dip \ rdr - — 7ia 2 , 
oo ^ 

where o = MIA is the area density of the island. The ratio Mil 
which appears in Equation 4 is thus equal to (867/608)a 2 = 
1.426a 2 . The equations of motion (Equation 3 and Equation 4) 
have been solved for k = 2,3,4 and increasing values of b until 
the island starts moving in the negative y direction and is lost. 
In Figure 3 the angle of motion 9 = dxoXdxsidYldX) is plotted as a 
function of the parameter b. The initial coordinate Yq of the tip 
along the slow scan direction was randomly chosen, with hardly 
any influence on the final results, except in the threshold region 
where the islands can be lost (and no points can be plotted). In 
all cases the direction of motion 6 initially decreases with 
increasing b and, again 9 — > 90° when b — > 0. However, the 
trend of the function 9(Z?) suddenly changes when b reaches a 
certain value (b = 0.5, 0.35 or 0.25 when k = 2, 3 or 4). In order 
to understand what happens at these points, we have also 
plotted the angular velocity of the particles, d<3)/dN, as a func- 
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Figure 3: Angle of motion 9 of 2/c-branched symmetric islands as a 
function of the distance b be-tween consecutive scan lines (in units of 
the length parameter a in the text), k = 2 (squares), k = 3 (circles) and 
k = 4 (triangles). 




tion of b (Figure 4). The critical values of the parameter b 
correspond to the onset of rotations of 1 807& angles per scan 
line. Beyond these critical values the angular velocity remains 
almost constant and the function 9(Z?) slightly increases 
(Figure 3). When k = 3 and 4 two further critical values of b are 
found (b = 0.9 and 0.35 respectively), corresponding to rota- 
tions of 2 x 1807A: angles per scan line. On the other hand, 
when b is small enough, the angular velocity d$>ldN becomes 
negligible: Rather than rotating, the islands simply 'wobble' 
like the nanorods. 

Discussion 

The predictions of the collisional model have been experimen- 
tally verified with gold nanospheres and nanorods manipulated 
on silicon oxide under ambient conditions by tapping AFM 
[6,7]. Furthermore, we have also observed that, at least in the 
case of the nanospheres, the model goes beyond the restrictive 
hypothesis that the particles are immediately stopped after being 
released by the tip. This has been shown by numerical simula- 
tions, where a 'mean free path' d of the nanoparticles was intro- 
duced. If the friction force between particle and substrate 
decreases, and consequently the distance d increases, then the 
pathway of the nanoparticle fluctuates more and more, but the 
form of the function 9(Z>) remains essentially unchanged [8]. 
Another important point is the following. In many commercial 
AFMs, the tip follows a zigzag scan path rather than a raster 
scan path. This leads to significant variations in the impact 
angles between the tip and particles and to a dependence of the 
direction of motion on the initial position of the particles along 
the fast scan direction x [7]. Nevertheless, at least in the case of 
nanospheres, one of the previous conclusions holds: The angle 
of motion 9 — > 90° when b — > 0 (in the case of a zigzag scan 
pattern, b can be taken as the distance between the starting 
points of parallel scan lines). Altogether, these observations 
suggest a general strategy for manipulating relatively large 
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nanoparticles, i.e., in the order of or larger than the tip radius. 
Provided that the density of scan lines is high enough, the direc- 
tion of motion of the particles can be tuned by orienting the fast 
scan direction x of the AFM perpendicularly to the desired 
direction of motion. This is much easier and more reliable than 
aligning the tip and moving it towards the COM of each 
nanoparticle, as is usually done. The rotational effects predicted 
by the collisional model have not yet been tested experimen- 
tally. A good benchmark would be the flower- shaped Sb islands 
first manipulated by Ritter et al. on HOPG and M0S2 [9]. 
Possible discrepancies between theory and experiment 
concerning the direction of motion and angular speed of the 
islands could be related to the friction forces between island and 
substrate and even used to estimate these forces in further 
developments of the collisional model. Since Sb islands can be 
manipulated and the corresponding friction forces can be 
measured also in contact mode [10], the applicability of the 
model could also be tested under these different impact condi- 
tions. Controlling the direction of motion of arbitrarily shaped 
nanoparticles is important for the guided formation of nano- 
structures. An interesting analogy is found with AFM nano- 
lithography. In a recent paper we have shown that the patterning 
of amorphous polymers can be 'tuned' by varying the scan path 
of an AFM tip which scratches the polymer surface while scan- 
ning [11]. Linear and 'travelling' circular ripples were formed 
using a raster or a circular scan path, respectively. In the same 
way, a desired configuration of nanoparticles could be obtained 
by a proper choice of the scan pattern. 

Conclusion 

In conclusion, we have shown that the direction of motion of 
nanoparticles can be controlled by AFM in a variety of 
significant cases. The key parameter is simply the density of 
scan lines in the scan path of the probe tip. Orienting the fast 
scan direction perpendicularly (and not parallel) to the desired 
direction of motion is an efficient way for manipulating the 
nanoparticles. With a proper choice of the scan pattern, it may 
be possible to reorganize an ensemble of randomly distributed 
nanoparticles in a well-defined arrangement. 
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